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LATTICE QCD WITH LIGHT TWISTED QUARKS: FIRST RESULTS 
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ON BEHALF OF THE EUROPEAN TWISTED MASS COLLABORATION (ETMC) 

I report on first results of the ETMC obtained simulating lattice QCD with two degenerate 
flavours of Wilson twisted mass fermions, for which physical observables are automatically 0(a) 
improved. Recent improvements of the HMC algorithm allow simulations of pseudoscalar masses 
below 300 McV on volumes ■ 2L with L > 2 fm and at values of the lattice spacing a < 0.1 fm. 
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1. Lattice QCD 

Simulations of lattice QCD are exact on the 
given lattice up to statistical errors, but 
many systematic uncertainties have to be 
controlled to have a reliable comparison be- 
tween lattice QCD computations and exper- 
imental measurements. 

• Renormalization and lattice artefacts: the 
choice of the lattice QCD action enters in 
the regularization of the theory. For this 
reason it is important to choose a lattice 
action that enjoys simple renormalization 
properties and good scaling behaviour. 

• Light quark masses: it is crucial to use 
algorithms that allow simulations at light 
enough quark masses in order to eventu- 
ally bridge the remaining gap between the 
simulation and the physical points, using 
chiral perturbation theory^ (^PT). 

• Finite size effects: the volume has to be 
large enough in order to avoid systematic 
uncertainties stemming from the finite size 
of the system that is simulated. 

The final goal is to simulate two light quarks 
and two non-degenerate heavier quarks keep- 
ing under control all these systematic uncer- 
tainties. Here I summarize first results for 
two light flavours obtained using a theoreti- 
cally well founded lattice QCD action and a 
highly improved algorithm. 



2. Wilson twisted mass QCD 

The lattice QCD action has a gauge Sg and 
a fermionic Sp part. The so called Wilson 
twisted mass (Wtm) action^ is 



where 

= 



[7m (V^ 



mo (2) 



is the Wilson-Dirac operator, toq is the un- 
twisted (bare) quark mass, fi the (bare) 
twisted quark mass and T3 the third Pauli 
matrix acting in flavour space. While parity- 
isospin is violated only at O(a^), the Wtm 
action has several main advantages: 

• The twisted mass provides a sharp infrared 
cutoff for low eigenvalues^. 

• Automatic 0(a) improvement when mg is 
properly tuned"^. 

• The renormalization pattern is in many 
cases greatly simplified^'^. 

Note that also with the standard Wilson 
{fi = 0) formulation 0(a) improvement can 
be obtained, but a whole set of parameters 
needs to be tuned to achieve this, while with 
Wtm only one parameter needs to be tuned, 
Too- To be more specific automatic 0(a) im- 
provement is at work if the value of mo is 
tuned to be the critical mass nir. In all the 
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results presented here the critical mass is ob- 
tained tuning the PCAC mass to zero and 
the role of the quark mass is played by the 
twisted mass ^. For recent summaries on 
Wtm see the reviews [6, 7] and references 
therein. 

2.1. Pseudoscalar decay constant 

The pseudoscalar decay constant /ps is an 
important quantity to compute for phe- 
nomenological applications (e.g. for the ex- 
traction of Vus from the ratio fx /fir), and 
can be used to fix the lattice spacing in 
physical units. With Wtm this computation 
presents a set of advantages: no improve- 
ment coefficients are needed (namely Cgw and 
ca) due to automatic 0(a) improvement and 
moreover no renormalization factors (Za) are 
needed^. All the uncertainties related to the 
estimates of these parameters are simply ab- 
sent, /ps and many other physical quan- 
tities have been studied'^ in the quenched 
model confirming that automatic 0(a) im- 
provement is at work and that the O(a^) ef- 
fects are small. 
3. Algorithms and phase 
structure 

To perform simulations at light quark masses 
with dynamical fermions, a substantial 
amount of computer resources is required. 
At the Lattice 2001 conference the situation 
was rather dramatic^. The low values of 
the quark masses, the large volume and the 
small lattice spacing reached in our work be- 
came possible only due to recent algorithmic 
developments^"'^^. In particular, for all the 
Nf = 2 results reported here we have used 
our variant of the HMC algorithm described 
in Ref. [11]. In continuum QCD with Nf > 1 
chiral symmetry is spontaneously broken at 
zero quark mass. The order parameter of this 
phase transition is the chiral condensate that 
is discontinous around the massless point. 
On the lattice Wilson-like actions break ex- 
plicitly chiral symmetry even in the mass- 
less limit. It is then a natural question, if 



one is interested in simulations of light quark 
masses, to understand the chiral phase struc- 
ture with Wilson-like quarks. Our collabora- 
tion has shown in a set of publications^^'^'^'^'* 
that there is a non-trivial phase structure at 
finite lattice spacing with a first order phase 
transition close to the chiral limit. These re- 
sults are in agreement with one of the pos- 
sible scenarios predicted^^'^^ by lattice-xPT. 
In this scenario the first order phase transi- 
tion line extends in the twisted axis direc- 
tion until a value of the twisted mass p, of 
order a^. As a consequence there is a mini- 
mal value of the quark mass that can be sim- 
ulated at fixed lattice spacing. Both the oc- 
curence of this scenario and the precise value 
of p, depend on the details of the lattice ac- 
tion, in particular on the choice of the gauge 
action Sq- The gauge actions so far studied 
by our collaboration can be parameterized by 

SG^f3[bo J2 (1- + 

x-.ij,<v 

+ &i(l- Vx2(.x;^,^.))] (3) 

with the normalization condition 6o = 1 ~ 
86i. pi^i indicates the standard plaquette 
and pi^2 a planar 1x2 rectangular loop. 
In particular it has been shown^^ that the 
value of fl decreases if bi is properly chosen. 
The choice of our collaboration was then to 
use bi = —1/12 which corresponds to the so 
called tree-level Symanzik (tlSym) improved 
gauge actiou"'^^. This choice compromises be- 
tween a smaller value for fl than obtained 
with 6i = and avoiding problems with big 
scaling violation possibly induced by a too 
large value of |6i|. 

4. First results 

4.1. Setting the scale 

In order to compare results obtained at dif- 
ferent lattice spacings and to translate the 
simulations results in physical units, it is nec- 
essary to set the scale, i.e. to give the lattice 
spacing in physical units. Many choices are 
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Table 1. The parameters for the simulation. In or- 
der to convert to physical units, I take the value of 
ro/a as computed at the minimal values of the twisted 
mass simulated, that can be read from Table 2. 



/3 


L3 X T 


ro/a(a/iniin) 


a [fm] 


3.9 


24^ X 48 


5.184(41) 


~ 0.096 


4.05 


32^ X 64 


6.525(101) 


~ 0.075 



Table 2. Current status of the simulations at /3 = 3.9 
and /3 = 4.05 using the parameters in Table 1: cur- 
rent estimates of the pseudoscalar and the renormal- 
ized quark masses, the number of trajectories gen- 
erated after 1500 trajectories of thermalization. 





fjR, [McV] 


mps [MeV] 




3.9 


~ 18 


~ 270 


5000 




~ 29 


~ 350 


5000 




~ 38 


~ 400 


5000 




~ 45 


~ 430 


5000 




~ 67 


~ 530 


5000 


4.05 


~ 18 


~ 270 


4200 




~ 34 


~ 370 


3000 



possible and some of them are under investi- 
gation by om' collaboration. Here I present 
results obtained by using as an hadronic in- 
put the force between two static quarks at 
a certain intermediate distance ro^^. While 
this quantity can be measured on the lat- 
tice very precisely, it has a rather uncertain 
phenomenological value. Nevertheless it can 
certainly be used to compare simulations at 
different lattice spacings. In this proceeding 
contribution to translate our results in phys- 
ical units I use the value vq = 0.5 fm. 

4.2. Pseudoscalar mass and 
decay constant 

The current status of our simulations are 
summarised in Tables 1 and 2. In Table 1 
the parameters of our simulations are pro- 
vided, and in Table 2 current estimates of 
the pseudoscalar mass mps and the renor- 
malized quark mass in physical units are 
given. A rough estimate of the missing renor- 
malization factor can be obtained from our 
data, using as inputs the available lattice 
estimates for the strange quark mass and 
'm-s/'m-u,d — 26. In Figure 1 the quark mass 
dependence of the pseudoscalar mass squared 
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is shown. The dependence is to a good ap- 
proximation linear, and the results at differ- 
ent lattice spacings do not show appreciable 
cutoff effects. I recall that the twisted quark 
mass in this plot still needs to be rcnormal- 
ized. In Figure 2 I show the dependence 
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Fig. 1. Quark mass dependence of the pseudoscalar 
mass squared. 

of the pseudoscalar decay constant on the 
squared pseudoscalar mass. The rather con- 
sistent values obtained at the two different 
lattice spacings is a very promising feature 
as far as the control on the systematic un- 
certainties of this physical quantity is con- 
cerned. For both the plots I have used xPT 
in order to correct for finite size effects^". In 
a forthcoming publication^^ we will analyze 
these data according to xPT. Here in Fig. 2 
I simply show for illustration a linear fit to 
the 5 lightest pseudoscalar masses, and the 
phenomenological estimate^'^^ of the pseu- 
doscalar decay constant in the chiral limit. 

5. Conclusions and outlooks 

The control of systematic uncertainties on 
lattice computations is a prerequisite in or- 
der to compare them with experimental mea- 
surements. I have summarized here some 
properties of Wtm, together with some first 
results with Nf = 2 light dynamical quarks. 
These results (see also the proceedings^^ at 
the Lattice 2006 conference) are rather en- 
couraging and suggest that the systematic 
uncertainties can be controlled simulating 
improved lattice actions, with the newly de- 
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Fig. 2. Dependence of the pseudoscalar decay con- 
stant on the squared pseudoscalar mass for 2 different 
lattice spacings. As a comparison is also shown the 
phenomenological value (a) in the chiral limit, and 
a linear fit to the 5 lightest pseudoscalar masses. 

veloped algorithms on currently available su- 
percomputers. 
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